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Experimental  and  CFD  Studies  of  Non-Thermal  Plasmas  for  NO,  Control 

M. J.  Nusca,  P. A.  Sagear*,  A.W.  Miziolek 


U.S.  Army  Research  Laboratory,  AMSRL-WM-PC 
Aberdeen  Proving  Ground,  MD  21005-5066  USA 

We  present  some  early  results  of  a  comprehensive  experimental  and  modeling  study  to  understand  the  driving 
physical  and  chemical  processes  involved  in  non-thermal  plasma  control  of  toxic  gases  such  as  the  NO,s. 


1.  Introduction 

Non-thermal  plasmas  (NTPs)  represent  a 
promising  technology  for  the  destruction  of  various 
toxic  gases  entrained  in  air  streams.  One 
diallenging  example  of  the  need  for  hazardous  gas 
control  is  the  jet  engine  test  cell  (JETC)  facility 
whCTe  vety  large  amounts  of  NO^  are  emitted  along 
with  various  hydrocarbon  gases  on  an  episodic 
basis.  The  high  flow  rates  as  well  as  the  wide 
dynamic  range  of  pollutant  concentrations  make 
this  a  particularly  difficult  problem  for  emissions 
control. 

In  order  to  investigate  the  ability  of  NTPs 
to  dissociate  NO^,  we  have  constructed  a  small 
scale  channel-flow  reactor  with  wall-mounted 
plasma  discharge  sites.  Toxic  gaseous  mixtures, 
including  NO^s,  will  be  entrained  in  an  airstream 
and  passed  through  the  reactor.  We  will  measure 
the  concentrations  of  various  radicals  produced  in 
the  reactor  by  the  destruction  of  NO,^.  This  data 
will  be  compared  to  results  generated  using  a 
computational  model  of  the  reactor.  Our  approach 
consists  of  closely  coupled  experimental  and 
computational  efforts  to  better  understand  the 
detailed  chemistry  and  physics  as  well  as  to 
simulate  and  scale-up  plasma  reactors. 

2.  Experimental 

The  laser  induced  fluorescence  (LIF) 
experiment  is  similar  to  that  used  by  other 
investigators  [1,2]  so  only  a  brief  outline  of  the 
instrument  is  given  here.  The  laser  system  for  both 
LIF  and  planar  LIF  (PLIF)  experiments  for  OH 
detection  consists  of  a  Lambda  Physik  Compex  1 02 
XeCl  excimer  laser,  a  Scanmate  2  dye  laser  using 
Coumarin  540  dye,  and  a  second  harmonic 
generator. 

For  single-point  LIF  measurements,  the 
laser  system  is  tuned  to  the  (1-0)  band  of  the 


-X^n  system  around  283  nm  and  focused  into  the 
plasma  reactor.  The  red-shifted  fluorescence  is 
centered  near  3 1 5  nm  (Fig.  1 )  and  is  collected  with 
a  two-lens  condenser  system  (Fig.  2).  The 
fluorescence  is  passed  through  a  broad-band 
interference  filter  centered  at  310nm  (50  nm 
bandwidth)  before  detection  by  a  PMT  (Phillips 
XP20I8B  or  Hamamatsu  R928).  After  passing 
through  a  preamplifier,  the  signal  is  averaged  with 


Figure  I .  Energy  Level  Diagram  for  OH  radical. 

a  SRS  SR250  boxcar  integrator  and  stored  in  a  PC. 
The  laser  intensity  is  measured  with  a  Scientech 
AC25HD  calorimeter  corihected  to  an  A/D 
converter  and  stored  in  the  PC.  Dispersed 
fluorescence  spectra  are  measured  by  focusing  the 
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Figure  2.  Schematic  Diagram  for  PLIF  studies. 


fluorescence  into  an  Acton  Research  Corporation 
S300i  monochromator  and  detecting  the  dispersed 
light  with  either  the  Hamamatsu  PMT  or  a 
Princeton  Instruments  IRY-1024G  intensified 
photodiode  array. 

PLIF  measurements  of  OH  are  made  by 
forming  a  laser  sheet  with  two  cylindrical  lenses, 
one  n^ative  and  one  positive  focal  length  (Fig.  2). 
The  sheet  is  passed  through  the  plasma  reactor  and, 
after  passing  through  an  interference  filter,  the  red- 
shifted  fluorescence  is  detected  with  an  intensified 
CCD  detector  (Princeton  Instruments  ITE/CCD- 
576G/RB-E)  equippefi  with  a  UV  Nikkor  camera 
lens.  The  transverse  intensity  of  the  laser  sheet  is 
simultaneously  monitored  by  using  a  fused  silica 
microscope  slide  to  direct  a  fraction  of  the  laser 
sheet  to  a  cell  containing  a  highly  concentrated  R6G 
dye  solution.  Another  camera  (Cohu  47 1 0)  detects 
an  image  of  the  resulting  fluorescence.  The  images 
are  stored  in  a  PC  for  off-line  image  processing 
later. 

The  plasma  cells  use  both  a  point-plane 
corona  discharge  and  a  silent  discharge  geometry 
and  were  constructed  at  LANL.  The  power  supply 
provides  high  voltage  pulses  up  to  50  kV  in 
amplitude  and  about  25  ns  in  duration.  Current  and 
voltage  profiles  are  simultaneously  measured  to 
monitor  the  energy  density  supplied  to  the  reactor. 

In  addition  to  the  detection  measurements 
on  OH  radicals  as  discussed  above,  temperature 
measurements  on  the  plasma  can  be  made  through 
analysis  of  the  relative  intensities  of  fluorescent 
bands  obtained  by  scanning  the  frequency  of  the 
dye  laser  which  has  a  bandwidth  of  0.15  cm’'.  An 
important  step  in  this  process  involves  the 


determination  of  the  ground-state  energy  levels 
from  the  relative  fluorescence  intoisities  of 
transitions  involving  known  rovibronic  levels.  The 
population  of  a  vibrational  level  in  the  ground 
electronic  state  is  simply  related  to  the  signal 
intensity  by  the  expression 

(1) 


where  ^  is  the  signal  intensity  and  /,  is  that  of  the 
laser.  Bjf  is  the  stimulated  absorption  coefficient  for 
the  transition,  and  is  the  fluorescence  yield  which 
can  be  calculated  from  stimulated  emission 
coefficients  and  nonradiative  decay  rates  due  to 
quenching  and  predissociation  [3,4].  The  Einstein 
coefficients  and  ground-state  rotational  energies  are 
all  readily  available  [5-8]. 

The  relative  populations  are  matched  with 
the  ground-state  energies  and  the  data  points  are 
plotted  and  fit  to  a  straight  line  according  to  Eq.  (2). 

N.=g.e 


The  transitions  of  interest  for  temperature 
measurements  are  in  the  (0,0)  band  of  the  A^S  - 
X^n  system  around  307  nm.  [3].  This  approach  is 
based  upon  the  assumption  of  rotational 
equilibrium  of  OH  in  the  non  thermal  plasma;  we 
are  continuing  to  investigate  the  validity  of  this 
hypothesis. 

Our  measurements  include  the  detection  of 
other  radicals.  In  particular  O,  H,  and  N  atoms  can 
be  detected  in  the  plasma.  The  detection  schemes 
for  these  atoms  invariably  involve  the  use  of  two- 
photon  spectroscopy  due  to  the  difficulty  of 
generating  vacuum  UV  radiation.  The  energy  level 
diagrams  for  these  species  are  shown  in  Fig.  3  and 
are  based  on  previously  reported  work  [9-13]. 

Data  from  the  experimental  reactor  is 
forthcoming  and  will  be  reported  in  future  papers. 
One  of  the  goals  of  the  present  study  is  the 
comparison  of  th^e  data  with  a  computational  fluid 
dynamics  (CFD)  model.  This  model  is  described  in 
the  next  section  and  the  comparison  between  data 
and  the  model  will  be  reported  in  future  papers. 


3.  CFD 

A  computational  fluid  dynamics  code,  used 
to  simulate  multi-component,  gaseous,  chemically 
reacting  flows  in  various  gasdynamic  applications, 
has  been  modified  for  numerical  simulation  of  the 
plasma  reactor  system.  The  ARL-NSRG2  code  was 
written  at  the  Army  Research  Laboratory  for 
solution  of  the  2D/axisymmetric  Navier-Stokes 
equations  including  species  diffusion,  thermal 
conduction,  viscosity,  and  nonequilibrium  chemical 
kinetics  [14].  An  explicit  fourth-order  Runge-Kutta 
algorithm  is  utilized  for  time-integration.  The 
convective  and  transport  terms  are  resolved  using 
central-differencing.  Numerical  stiffness  due  to 
chemical  source  terms  is  mitigated  by  convergence 
acceleration  techniques.  Flow  turbulence, 
introduced  as  an  in-flow  condition,  is  modeled 
using  a  mixing  length  approach. 

The  electrons  added  by  the  plasma  through 
microdischarges  in  the  flow  field  are  modeled  using 
source  terms  in  the  species  mass  conservation 
equations.  These  discharges  promote  chemical 
reactions  in  regions  of  the  flow  field  interacting 
with  the  plasma.  An  initial  NO„  chemical  kinetics 
mechanism,  incorporated  into  the  code,  includes  9 
reactions  and  12  species.  The  reactions  considered 
consist  of  three  dissociation  reactions  (Oj,  N  j,  and 
HjO  to  0,N,  and  OH),  three  ozone-producing 
reactions,  and  three  reactions  that  produce  NO  and 
HNO3  from  0,N,0H,  and  NO2. 

Figure  4  shows  the  specific  configuration 
used  in  this  initial  study  to  test  the  CFD  code.  The 
inflow  is  uniform,  laminar,  and  consists  of  wet  air 
and  NOj.  Although  the  reactor  has  sides  (i.e.,  10 
cm  apart),  the  present  CFD  simulation  was  run  as 
2D  along  the  centerplane  of  the  reactor,  as  shown  in 
the  lower  part  of  Figure  4.  Plasma  discharge 
locations  were  placed  10  cm  apart  and  were 
assumed  to  discharge  at  a  frequency  of  1  kHz.  Both 
the  inflow  and  discharge  temperatures  were  300K. 
The  lower  part  of  Figure  4  shows  the  computational 
mesh,  100  cells  in  the  flow  direction  and  40  cells  in 
the  transverse  direction.  The  mesh  is  clustered  near 
the  reactor  walls  (i.e.,  plates)  in  order  to  increase 
accuracy  of  the  computed  wall  boundary  layers. 
The  mesh  is  uniform  in  the  flow  direction.  The  four 
dark  lines  are  not  physical,  but  are  added  to  indicate 
the  locations  of  the  discharges. 

Figure  5  shows  a  representative  result  of 
CFD  simulation.  The  top  part  of  the  figure  displays 


the  fully-developed  flow  in  the  channel  (reactor)  as 
velocity  vectors  overlaid  with  grey-scale  contours  of 
non-dimaisional  velocity.  Boundary  layers  develop 
on  the  reactor  walls  (i.e.,  small  vectors)  and  the  core 
velocity  is  about  0.9.  The  bottom  view  of  the  figure 
displays  grey-scale  contours  of  electron  mass 
fraction  at  the  instant  the  discharges  occur;  at  four 
locations  the  mass  fraction  of  electrons  is  unity 
(dark)  while  elsewhere  the  mass  fraction  is  null 
(white). 

Figure  6  shows  grey-scale  contours  of 
NO2  mass  fraction  (note  the  restricted  scale  of  mass 
fraction)  at  near  the  instant  the  discharges  occur. 
Depletion  of  NOj  occurs  at  the  discharge  sites  as 
NO2  reacts  with  O,  N,  and  OH  (formed  by  O  2,  N  2 
and  H2O  reactions  with  electrons)  to  form  NO  and 
HNO3.  The  bottom  view  of  the  figure  shows  the 
distribution  of  HNO3,  as  number  density,  along  the 
center  of  the  reactor  (i.e.,  x=0  to  50  cm,  y=.l  5  cm) 
at  various  times.  The  discharge  is  occurring  every 
0.001s.  The  slow  buildup  of  HNO3  is  indicative  of 
reactions  occurring  to  breakdown  NO2  throughout 
the  reactor. 

The  major  assumptions  used  in  the  present 
CFD  simulations  are  laminar  flow,  certain  reduced 
kinetics  mechanism  and  rates,  and  certain 


(n^2) 


(n=1) 


205  nm 


(n=0) 


"0 


656  nm 


211  nm 


669  nm 


*P 


226  nm 


845  nm 


H 


'S 

N 


jp 

O 


Figure  3.  Energy  level  diagrams  for  H,  N,  and  O 
atoms. 


assumptions  about  diffusion  and  thermal  properties 
of  the  species.  The  results  shown  in  Figures  5  and 
6  are  influenced  by  these  assumptions.  The  CFD 
code  has  been  setup  for  the  reactor  problem  and  is 
now  available  to  test  the  dependence  of  the  results 
on  the  assumptions  and  do  some  further  work.  The 
introduction  of  turbulence  in  the  inflow,  at  the  walls 
and  with  flow  obstructions  will  cotainly  effect 
these  results. 
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DEMONSTRATION  PROBLEM  SETUP 


PLASMA  DISCHARGE  LOCATIONS  (E=25  MJ/cm^,  (o=1  kHz,  T=300K) 
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Sample  Chemistry  Model,  12  species:  e",  Oj,  Nj,  O,  N,  0„  NO,  NO,,  H,  OH,  H,0,  HNO, 

9  reactions:  0/N,/H,0  dissociation,  ozone-creation,  NO,  reduction. 
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Figure  4.  Probiem  Setup  and  Computational  Mesh. 
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Figure  5.  Computed  Flowfields  with  and  without  Plasma  Discharges. 


Figure  6.  Computed  NO^  Flowfields  and  HNO,  Distributions. 


